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Abstract 

PAMELA and ATIC recently reported an excess in e ± cosmic rays. 
We show that if it is due to Dark Matter annihilations, the associated 
gamma-ray flux and the synchrotron emission produced by e ± in the 
galactic magnetic field violate HESS and radio observations of the 
galactic center and HESS observations of dwarf Spheroidals, unless 
the DM density profile is significantly less steep than the benchmark 
NFW and Einasto profiles. 
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1 Introduction 



Evidence for the existence of Dark Matter comes from a number of astrophysical and cosmolog- 
ical probes [1], yet no direct uncontroversial detection has been made so far, and its properties 
remain largely unknown. The excesses in cosmic ray e spectra suggested by the recent ob- 
servations of PAMELA and ATIC-2 in the energy range 10 GeV - 1 TeV can be tentatively 
interpreted as due to DM annihilations in the galactic halo. Indeed, the PAMELA satellite 
reported [2] an excess in the positron fraction (e + /(e + + e~)), with respect to the expected as- 
trophysical background above 10 GeV. The spectrum features a steep rise up to 100 GeV, the 
highest energy currently probed by the experiment. No excess is instead seen in the antiproton 
flux [3]. The balloon-borne experiment ATIC-2 [4] reported the detection of a peak in the 
undistinguished flux of positron and electrons, in the energy range between 300 and 800 GeV. 
The two signals are compatible and, at low energy, the precise data from PAMELA confirm 
previous hints from other experiments [5]. 

Of course caution should be used when interpreting the data, as the background flux from 
conventional astrophysical processes (e.g. spallation of cosmic rays on the interstellar gas) carry 
big uncertainties, as large as one order of magnitude, so that even the very evidence for an 
excess in PAMELA data could possibly be jeopardized [7]. A more likely possibility is that 
the observed effects are due to some single astrophysical object, such as a pulsar (a foreground 
signal for astrophysics, but an annoying distraction for Dark Matter searches), or a collection 
of them. Indeed it is expected that pulsars produce a power law spectrum of mostly electron- 
positron pairs, with a cut-off in the multi-TeV range. The known nearby pulsars Geminga 
and B0656+14 are the main candidates, but unknown and past pulsars can contribute to the 
integrated flux [8]. Work is underway to assess more precisely the features of the expected 
fluxes from pulsars and discriminate from the DM hypothesis [9]. 

If however Dark Matter annihilations are at the origin of the observed anomalies in cosmic 
rays, than the data allow to determine the required properties of the Dark Matter: mass M, 
annihilation cross section av and main annihilation mode. The analysis in ref. [10] studied 
systematically which of these can fit the data, finding three main classes of models. 

a) on the basis of the e + and p data from PAMELA, the Dark Matter can be: 

al) a particle that dominantly annihilates into leptons, with no strong preference for the 
mass, if above a few hundred GeV; 

a2) a particle that annihilates into W, Z or higgses and that has a mass > 10 TeV. 
An example in this class is the fermion quintuplet in the predictive Minimal Dark 
Matter model [6], where DM has mass M = 9.6 TeV and annihilates into transversely 
polarized W + W~: this model predicted the PAMELA excess and predicts that the 
ATIC peak is not there; 

b) adding the peak from ATIC-2, a clear indication for the mass emerges: DM has to be a 
particle with mass ~ 1 TeV that dominantly annihilates into leptons. We will exemplify 
this class of models referring to a generic candidate with M = 1 TeV and annihilations 
into p, + \i~ . 

The upcoming results of ATIC-4 [11], PAMELA, or the first data from the Fermi LAT 
calorimeter [12] or Air Cherenkov Telescopes [13] can soon check if a peak is really present in 
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the e + + e~ spectrum just below 1 TeV: if the peak is there b) is favored and a) is excluded; 
if instead the peak is not there, then a) is favored and b) excluded. Models with M< 1 TeV 
appear to be already disfavored. 

For what concerns the magnitude of the annihilation cross section, the large flux above the 
background in the PAMELA and ATIC data indicates a very large av (see fig. 9 of [10]). For 
instance, for a candidate in class b), a value of the order of few 10~ 23 cm 3 /sec is needed to fit 
the data. This is much larger than the typical thermal cross section av = 3 • 10~ 26 cm 3 /sec 
suggested by the cosmological DM abundance. As discussed in [10, 22], the two values can be 
reconciled if a Sommerfeld enhancement is at work: this effect, in fact, depends on the DM 
velocity in an important way and so it would be present at v ~ 10 -3 (the typical velocity of 
DM particles annihilating in the galactic halo at the present time) and reduced or absent at 
v ~ 0.2 (the velocity at decoupling). More precisely, the enhancement of a non-relativistic 
s-wave DM annihilation [19, 20, 21, 22] can be approximatively characterized in terms of two 
critical velocities u mm and t> max as follows: 



In terms of particle-physics parameters, assuming that the long-range force that gives rise to 
the Sommerfeld enhancement is a vector with mass My <C M and gauge coupling gy to DM, 
one has t> max ~ <7y/4. The value of t> min is t> min ~ M v /M unless a (DM DM) bound state 
with small binding energy E B is present; in such a case the Sommerfeld effect grows down to a 



smaller f min « J Eb/M. In the exemplar model of class al), the enhancement is automatically 
present (via the exchange of weak gauge bosons). Extra states are instead required in class 
b), so that [10] suggested that DM might be charged under an extra U(l), proposing a specific 
model. More proposals along these or different lines have followed [22, 14]. Alternatively one 
can invoke either non-thermal DM or very large boost factors [23]. However, the latter can 
only arise in rather exotic scenarios (e.g. DM mini-spikes around black holes [17, 18]), but not 
in the framework of DM subhalos with realistic properties [16]. 

Given these tantalizing but surprising hints of Dark Matter annihilations in the charged 
particle signals, it is now crucial to consider the constraints on this interpretation that come 
from the photon fluxes that necessarily accompany such charged particles. These photon fluxes 
are produced: 

i) directly as a product of the DM annihilations themselves (mainly from the bremsstrahlung 
of charged particles and the fragmentation of hadrons, e.g. tt°, produced in the annihila- 
tions), at energies comparable to the DM mass M, i.e. in the 7-ray energy range of tens 
of GeV to multi-TeV. 

ii) at much lower energies, e.g. radio to visible frequency, by the synchrotron radiation emit- 
ted in the galactic magnetic field by the electrons and positrons produced by DM anni- 



The best targets to search for these annihilation signals are regions with high DM densities, 
such as the Milky Way Galactic Center (GC), the Milky Way Galactic Ridge (GR) and the 
Sagittarius Dwarf spheroidal satellite galaxy (Sgr dSph). The predicted photon fluxes can then 
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Figure 1: Shape of DM density (left) and magnetic field (right) profiles discussed in the text, 
as a function of the galactocentric coordinate r. 



be compared with observational data, in order to rule out combinations of astrophysical and 
particle physics parameters that violate observational constraints. 

The aim of this paper is to compare the regions suggested by the PAMELA (and ATIC) 
data in the plane of annihilation cross section and DM mass (av, M) with those excluded by 
photon observations. We perform the analysis for arbitrary values of M and for several different 
primary annihilation modes. We take into account different choices for the main astrophysical 
unknown ingredients: the galactic DM density profiles and the galactic magnetic field. In 
section 2 we discuss bounds from gamma-ray observations, mainly performed by the HESS 
experiment. Section 3 discusses bounds from lower energy photons radiated by the e ± . 

2 7 ray observations 

We start by considering the 7-ray fluxes produced by DM annihilations directly. Since DM is 
neutral, a tree-level annihilation into 7's is of course not possible, thus the flux is the sum of 
various effects that arise at higher order in a em : i) a continuum at lower energies produced 
by the bremsstrahlung of charged particles and the fragmentation of hadrons produced in the 
annihilations; ii) a line at E m M produced by one-loop effects; iii) possibly a continuum at 
E just below M produced by three-body annihilations [15]. Infrared divergences in the total 
annihilation rate cancel among i) and one loop corrections without photons in the final state, 
and these contributions are separately gauge invariant in the energy ranges where they are 
separately relevant. The details of contributions ii) and iii) are model dependent, so that we 
only consider the contribution i). 

The differential flux of photons from a given angular direction dQ is 

rf$ 7 = lr p% dNf 1 chfp(r)\ 2 

dVtdE 2 4ttM£ m y K /f dE ' ./line— of —sight \ Pq / 

where r « 8.5 kpc is the distance of the Sun from the galactic center, p = 0.3 GeV/cm 3 is the 
DM density at the location of the solar system and / runs over all the 7-ray producing channels 
with annihilation cross section (crv)f and individual spectrum dN^/dE. The adimensional 
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MW halo model 


r s in kpc 


p s in GeV/cm 3 


J (io- 5 ) 


NFW [24] 


20 


0.26 


15 • 10 3 


Einasto [26] 


20 


0.06 


7.6 ■ 10 3 


Isothermal [27] 


5 


1.16 


13 



Table 1: Parameters of the density profiles for the Milky Way discussed in the text and cor- 
responding value of J for AQ = 10~ 5 . In all cases we imposed the normalization p(r ) = 
0.3 GeV/cm 3 . 



quantity J encodes the astrophysical uncertainty. When observing a region with total angular 
size AQ the factor J dQ gets replaced by J ■ AQ = f AC J dQ. 



In order to compute the flux, one thus has to specify the DM density profile as a function 
of the galactocentric coordinate, p(r). It is commonly assumed that DM follows a 'universal' 
Navarro, Frenk and White (NFW) profile [24] 

PNFWO) = 7"^ T2 • ( 3 ) 



+ 



It is however unclear whether this analytic formula, obtained by fitting DM halos in iV-body 
simulations, actually captures the behaviour of the density profile down to the innermost re- 
gions. In fact, although there is general agreement on the shape of profiles at large scales, 
profiles steeper than NFW, with p(r) oc r" 12 at radii much smaller than the virial radius, have 
been found to provide a better fit to simulated halos [25]. This claim has been subsequently 
challenged, e.g. by Ref. [26], where it was found that the so-called Einasto profile 



PEinasto^J = Ps ' exp 



a \\r 



a = 0.17 (4) 



should be preferred, since the profiles of simulated halos appeared to become shallower and 
shallower towards the Galactic center, without converging to a definite power-law. 
Finally, a truncated isothermal profile 

Piso(r) = (5) 



is sometimes adopted as a benchmark, since it is representative of 'shallow' DM profiles [27]. 

In Table 1, we show the parameters of the aforementioned density profiles (plotted in fig. la) 
for the case of the Milky Way and the value of the quantity J averaged over a solid angle 
AQ = 10~ 5 str, corresponding to the angular resolution of gamma-ray experiments such as 
HESS and Fermi LAT. Note that, for the Einasto profile, we choose a value of r s = 20 kpc 
representative of the results of Ref. [26] for different simulations; the value of J (AQ = 10~ 5 ) is 
only a factor ~ 2 smaller than the one for NFW. 

We stress that aside from these uncertainties, the DM distribution is further complicated 
by a number of physical processes that are not accounted for in most numerical simulations, 
such as the presence of a supermassive black hole that dominates the gravitational potential 
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within 1 pc from the Galactic center, and a stellar cusp that inevitably interacts with the DM 
fluid, a circumstance that makes it difficult to accurately estimate the DM profile in the central 
region (see e.g. Ref. [28] and references therein). We do not include these model-dependent 
processes in the following. More generally, the extrapolation of the numerical Dark Matter 
profiles in the regions very close to the Galactic Center ( < 10 pc) is of course to be taken 
with simulations cannot resolve small radii. As we will see, however, some of the 

constraints come from regions as large as (9(100) pc (the size of the Galactic Ridge region, for 
instance, or of the Sagittarius Dwarf galaxy) where the impact of these uncertainties is much 
less important. Keeping in mind these remarks and possible caveats, in the following we will 
discuss the astrophysical constraints for different choices of the DM profile. 

2.1 7- ray observations of the Galactic Center 

HESS observations in the direction of the Galactic Center have revealed a source of Very High 
Energy 7-ray emission (HESS J1745-290) lying within 7" ± 14g tat ± 28g yst from the supermassive 
black hole Sgr A*, and compatible with a point source of size less then 1.2' [29]. The correspond- 
ing energy spectrum, shown in figure 2a, is well fitted by a power law dQ^/dE oc E~ 2,25±0M , 
over two decades in energy, and it has been confirmed by the MAGIC collaboration [31]. The 
EGRET experiment had actually previously reported the detection of a point source (3EG 
J1746-2851) within 0.2 degrees from Sgr A* [30]. However, a re-analysis based on photons with 
energies above 1 GeV has shown that the source is slightly offset with respect to the galactic 
center [32]. 

The possibility to interpret both sets of 7 observations (separately or at the same time) 
in terms of DM annihilations has been discussed e.g. in Refs. [33, 34, 35, 36, 37, 38, 39], and 
Ref. [40] discussed the prospects for the detection of DM with Fermi LAT. Here, we take a 
conservative approach and consider the observed gamma-ray emission as an upper limit to 
the DM annihilation flux, in order to test the compatibility with a DM interpretation of the 
PAMELA data. We compute the constraints in the av versus mass plane, by requiring that 
DM annihilation flux does not exceed (at 3a, in terms of the error bars quoted by the HESS 
collaboration) the observed emission at any data point. 

As an example of our confrontation with data, the left panel of fig. 2 shows the gamma-ray 
flux from the Galactic center (assuming a NFW profile) produced by the annihilations of 10 
TeV DM particles into W + W~ , the aforementioned sample model al), with an annihilation 
cross section av = 10~ 23 cm 3 /sec. This mimics the Minimal Dark Matter theory [21]. As one 
can see, the DM 7 flux does not exceed any of the HESS data points, and correspondingly this 
point in parameter space will lie in the allowed region of fig. 4 (see below). The left panel of 
fig. 2 also shows the superposition of the DM signal with a sample power-law background: while 
in this case one would conclude that the model is excluded because the summed flux exceeds 
the HESS observations by more than 3a at several data points, this conclusion would not be 
solid. Indeed, choosing a different background (e.g. lower in normalization) within its large 
uncertainties can re-allow the model. Adopting the criterion of comparing each single data 
point with the DM-only flux allows us to have more conservative and robust results, in the 
sense that they can only be made tighter by specific choices for the background, that we do not 
perform in order not to exclude more models than it is justified by current data. Moreover, we 
recall that the signal from DM annihilations that we consider here does not include all the model 
dependent contributions (see the discussion at the beginning of sec. 2). A full computation of 
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Figure 2: HESS observations of the Galactic Center (left) and Galactic Ridge (right) compared 
with the annihilation signals from our two sample models al) (left) and b) (right), assuming a 
NFW profile and an annihilation cross section of av = 1CT 23 cm 3 /sec, and with the sum of the 
annihilation signal and a possible astrophysical background flux. 



the gamma-ray spectrum and of its uncertainties to be inserted in a global fit is possible only 
in specific theories. For instance for the case of Minimal Dark Matter it has been obtained in 
Ref. [21]. These contributions can change somewhat the shape of the DM signal and bring it 
closer to the shape of the observed spectrum, even in absence of a power-law background. For 
this additional reason, it is apparent that it would be wrong to exclude a model such as the 
one illustrated in the left panel of fig. 2. 

The HESS collaboration has also recently discovered a diffuse gamma-ray emission, cor- 
related spatially with the Galactic Ridge (GR), a complex of giant molecular clouds in the 
central 200 pc of the Milky Way [41]. Once point sources, including HESS J1745-290, are 
subtracted, the reconstructed gamma-ray spectrum for the region with galactic longitude 
—0.8° < i < 0.8° and latitude \b\ < 0.3° is well described by a power law with photon in- 
dex T = 2.29 ± 0.07 sta t ± 0.20 syst . In this region, the predicted DM signal is smaller than in 
a small cone pointing towards the Galactic center, but the astrophysical background is also 
significantly reduced, and the constraints are less sensitive to the slope of the DM density pro- 
file. Fig. 2b shows the HESS data and the signal in our sample model b), i.e. a candidate with 
M = 1 TeV annihilating to assuming a NFW profile. The cross section has been taken 

here to be av = 10~ 23 cm 3 /sec as in the left panel. The same discussion as above applies: 
what is seen here is that the DM 7 signal exceeds a couple of data points by more than 3a and 
therefore the model will lie within the excluded region. 

In figures 3 and 4 we show the results of the analysis of the data described above. The 
continuous blue lines shows our conservative bounds on the annihilation cross section av from 
HESS observations of the Galactic Center, and the dot-dashed blue lines show the comparable 
bounds from Galactic Ridge observations. Figs. 3 refer to DM annihilation into leptons, while 
fig.s 4 show the more 'traditional' DM annihilation modes into W + W~ , bb and ti. Barring the 
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Sgr dSph halo 


Parameters 


Core/Scale radius 


J (2- 


io- 5 ) 


J {2- IO" 5 ) 


Small core [42] 


v a = 13.4 km/s 


r c = 1.5 pc 


31 • 


10 3 


74 10 24 GeV 2 /cm 5 
2.46 10 24 GeV 2 /cm 5 
0.32 10 24 GeV 2 /cm 5 


NFW [50] 


p s = 5.2 GeV/cm 3 


r s = 0.62 kpc 


1 • 


10 3 


Large core [44] 


v a = 22.9 km/s 


r c = 0.23 kpc 


0.14 


■ 10 3 



Table 2: Parameters of the density profiles for the Sagittarius Dwarf galaxy discussed in the 
text and the corresponding value of J(Afl) (normalized by convention in terms of the solar 
quantities r and p Q , as in eq. (2)) for AQ = 2 • 10~ 5 . For reference, the value of the rescaled 
J~(AQ) = TqPq J(AQ) is also given. 



possibility of boost factors or Sommerfeld enhancements different for 7 and e observations, 
we see in fig.s 3 that the green regions that can fit the PAMELA anomaly (and the red regions 
that can also fit the ATIC anomaly) are excluded for masses M > 300 GeV, by two orders of 
magnitude if DM follows the NFW density profile, by an order unity factor if DM follows the 
Einasto profile, and are allowed if DM follows the isothermal profile (somewhat disfavored, 
however, by iV-body simulations). 

In fig.s 4 a similar situation holds. The green PAMELA bands are here truncated because low 
DM masses do not allow a good fit to the anti-proton data; the truncation is conservatively 
put at 1 TeV, but masses up to multi-TeV still do not give a good fit, as discussed in the 
Introduction (see [10] for the full analysis). 



2.2 Sagittarius dwarf spheroidal galaxy 

Dwarf spheroidal galaxies are among the most DM-dominated structures, so that they allow 
to search for 7 ray signals of DM annihilations with minimal astrophysical backgrounds. In 
particular, HESS has observed the Sagittarius Dwarf galaxy [42], a satellite of the Milky Way 
which is located at a distance of d = 24 kpc from the Sun. The satellite is thought to be in the 
process of being disrupted by multiple passages through the Milky Way disk, and the fact that 
it still exists is taken as an indication of the existence of a substantial amount of Dark Matter 
in it. 

The DM density profile in Dwarf Galaxies is uncertain as much as the one in the Milky 
Way, with which it might have some correlations. For Sgr dSph we consider the possibilities of 
a cusped NFW profile [42, 44] with density given by eq. (3) and of the class of cored profiles 

P (r)= Vl (6) 

The normalization factors and the characteristic radii are reported in Table 2, where also the 
corresponding values of J, defined according to eq. (2), are given. The area of observation 
corresponds to an aperture angle of 0.14° i.e. to a size of AQ = 2 ■ 10~ 5 [42]. 

HESS has observed Sagittarius Dwarf for T ^ s = 11 h finding no 7-ray excess: the integrated 
photon flux is iV 7 < 85 at 3cr (iV 7 < 56 at 95% CL [42]). Hence an upper bound can be imposed 
on the annihilation cross section 

8tt M 2 AL 
av< — ~ Is y A n ~ rfAT (7) 



r obs r Q p Q JAQ J dE A ^ E ^ 



dE 
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where the effective area of HESS for observations at ~ 20° (Sgr dSph is located at 14° galactic 
latitude) A cS (E) ~ 10 5 m 2 in the range E>70 GeV is taken from [43]. 

The resulting bounds on av are shown as dashed blue lines in figures 3 and 4. The top 
rows of the figures assume a NFW DM density profile in Sgr dSph: the bounds are overall 
comparable or slightly less powerful than the bounds from the Galactic Center and Ridge. In 
all the lower rows we use for Sgr dSph a 'large core' profile, which gives the minimum 7 flux 
among the profiles considered in the literature. The bound becomes the most constraining 
one when the Milky Way profile is taken to be isothermal. We have not explored whether even 
smoother profiles of Sgr dSph can be designed (compatibly with observations) that can lift such 
bound. It is interesting to note that the typical velocity dispersion of DM in Dwarf Spheroidal 
galaxies is about lOkm/s [44], smaller than in our galaxy: thereby their constraint becomes 
stronger and dominant for models where light particles give a Sommerfeld enhancement down 
to a small v m i Q [22] . 

Notice that the regions suggested by PAMELA for light DM mass are not probed by HESS 
observations due to its high energy threshold. DM that annihilates into leptons tends to give 
most of the signal at 7 energies just below the DM mass: the forthcoming Fermi/GLAST 
7 observations are not expected to be very significant for our purposes, as they will extend 
the HESS observations down to lower energies but will overlap with the HESS observation at 
the > 100 GeV gamma energies suggested by the PAMELA/ATIC excesses. We now turn to 
radio-wave observations. 



3 Radio observations of the Galactic Center 

The e produced by DM annihilations within the galactic magnetic field radiate synchrotron 
radiation. The Galactic Center is presumably the best region to search for this effect, because 
of the large local value of the DM density and magnetic fields [45, 39, 35, 46, 47, 48, 38]. We 
first detail the necessary astrophysical and particle physics ingredients, and then move to the 
comparison with observations. 

The GC region contains a black hole with mass Mbh ~ 4.3 x 1O 6 M (see e.g. the recent 
Ref. [49] and references therein). This implies two length-scales: the Schwarzschild radius 
Rbh = 2GnM B h f& 4 x 10 _7 pc and the radius of the accretion region, i? acc = 0.04 pc, defined 
to be the region where the velocity flow due to the gravity of the black hole, v = — \J RBn/ r is 
larger than the random galactic motion, v ~ 10~ 3 . 

Assuming a constant accretion of the BH mass, Mbh ~ 5 • 10~ 12 M /sec, the matter density 
is given by p(r) = Mbh/^^v (r) oc r -3 / 2 . Assuming equipartition of the matter kinetic 
energy with the magnetic pressure, pv 2 /2 = B 2 /2, the magnetic field is B(r < R aC c) = V pv 2 = 
7.2 mG • (-Race/ r) 5//4 . Outside the accretion region, assuming that the magnetic flux is conserved, 
B{r > Race) scales as 1/r 2 , down to the typical constant galactic value B ~ /iG reached at 
r ~ 100i? acc . This defines the 'equipartition' magnetic field, plotted in fig. lb. Another 
possibility is that B stays constant inside the accretion region: this defines the 'constant' 
magnetic field, again plotted in fig. lb. As we will see, bounds from observations at lowest 
frequencies are robust, and mildly vary even when the magnetic field is varied within these 
extremal possibilities. 
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Next, we need to compute the number density n e (r,p,t) of the e generated by DM an- 
nihilations. We assume stationary conditions, spherical symmetry, and, in view of the large 
magnetic fields, we neglect diffusion and assume that synchrotron radiation dominates energy 
losses. Writing the injection term (numerical coefficients are given in particle-physics natural 
units in the following) Q = av(p 2 /2M 2 )(dN e ±/dE) as Q = 47rp 2 q and the number density as 
n e = 4Trp 2 f, f(r,p) obeys the equation 



d f . 1 9 r- 2,1 , • df 



v~ + 



PsynP f +Padv^- = Q (8) 



with energy losses 



or p A op L J op 



p d(r 2 v) . e 4 B 2 Ep 
3r 2 or 9nmj 



We assume p ^> m e , such that E ~ p and E syn ~ p syn . If one can neglect advection (because 
Padv <S P syn : this happens at r > _R acc and at large p), the above equation is solved as 

1 r°° o 2 N (E) 

n e (r,E)^-^ dE> Q(E>,r)=av-P-^- (10) 



where N e (E) is the number of e + or e generated with energy larger than E in one DM 
annihilation. 

Finally, we can now compute the synchrotron power W syn generated by the n e electrons and 
positrons in the turbulent magnetic field B 

dW syn _ V3e 3 B rn/ v s f°° T , 8tt 



* -6. m /0 n*) = *f JMOde-j^-i/i!) (») 

where 



^S^MH^) . (12) 



Anml G 

Reducing the magnetic field S, the spectrum of synchrotron radiation moves to lower energies, 
but the total energy into synchrotron radiation remains constant, until B becomes so small 
that other energy-loss mechanisms start to dominate. 
Inserting eq. (10) in (11) we find 

where the integral extends over the observed volume and p is obtained from eq.s (11) and (12) as 
p = ^Anm^v/B = 0.43 GeV(z//GHz) 1 / 2 (i?/mG)~ 1 / 2 . A lower B leads to a higher synchrotron 
flux at the low frequency we consider. 

We now move to the comparison with observations. Since the observed GC microwave spec- 
trum is harder than what DM decays can produce, the dominant bound is obtained considering 
the observation available at the lowest observed frequency, v = 0.408 GHz, performed by [51] 
in a region with full width half maximum of 4". The observation found an upper limit to the 
measured flux S = [y dW syn / dv) / (Airr 2 ^) < 2 10 -16 erg/cm 2 sec, that constraints from above the 
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flux in eq.13. The resulting bounds are plotted in fig.s 3 and 4 as red lines. What is seen is 
that this constraint excludes a large portion of the parameter space for NFW and Einasto DM 
profiles. The constraint extends to low DM masses (where the 7-ray bounds from HESS are not 
effective). The variation of the magnetic field between 'equipartion' and 'constant' in the inner 
region at r < R acc negligibly affects the bound, because the radio emission is predominantly 
produced by the outer region. We have also verified that for the relatively shallow profiles 
under consideration, the synchrotron self-absorption is negligible. 

The subdominant bound (purple lines) comes from the VLT observation at the larger in- 
frared/visible frequency, v = 0.5 10 5 GHz: S < 3 1CT 12 erg/cm 2 sec from a region with angular 
size 0.04" i.e. r < 0.0016 pc. It somewhat depends on the magnetic field profile, and it be- 
comes numerically significant only for spiked DM density profiles [38]. Similarly, observations 
at higher frequencies give possibly strong but not robust bounds [38] , that also strongly depend 
on the possibility of having an intense 'equipartition' magnetic field close to the Milky Way 
black hole. 

Finally, we need to consider the possible effect of advection, that we neglected so far. We 
numerically studied it, and, as it depends on the density profile, magnetic field, DM mass 
and annihilation mode we try to qualitatively summarize some general lessons rather than 
presenting a large number of specific plots. At VLT-like (and higher) frequencies, v ~ 10 14 Hz, 
synchrotron radiation is dominantly generated by e with energy above 1 GeV. The density 
of such higher energy e is negligibly affected by advection. Presumably, advection can be 
neglected also for the observation by [51] at the lower frequency v = 0.408 GHz, because the 
angular aperture of 4" corresponds to a region with size i? Q b s = 0.14 pc, which is larger than 
the (presumed) accretion radius R acc ~ 0.04 pc, inside which advection is a significant effect 2 . 
Advection can be neglected because the total luminosity is dominated by radii comparable to 
i?obs, unless one considers DM density profiles that grow at r — > more strongly than the NFW 
profile, leading to very strong constraints possibly affected by advection (and by the reduction 
of the Sommerfeld enhancement, as the DM velocity grows, becoming relativistic around the 
Schwarzschild radius). The Davies bounds can possibly be weakened if the accretion region is 
larger. 

4 Conclusions 

We explored the compatibility of the interpretation in terms of DM annihilations of the ex- 
cesses in the CR e ± spectra claimed by PAMELA and ATIC with observations of photons at 
gamma and radio frequencies, inevitably produced by brehmstrahlung (in DM annihilations) 
and synchrotron radiation (of e in galactic magnetic fields) respectively. 

Figures 3 and 4 summarize our results. The green (red) bands show the region allowed by 
PAMELA (PAMELA and ATIC combined), varying the e propagation models between the 
so-called min/med/max sets of parameters (see e.g. [21, 53] and references therein). We do not 
show the subleading < ± 20% experimental and background uncertainty. The regions shaded 
in blue (red) are excluded by our fit of gamma (radio) observations. One sees that the two 
data-sets are incompatible if the DM density profile is Einasto, NFW or steeper: the gamma 

2 Advection typically reduces the e ± density just below i? acc and e ± accumulate in the inner region. 
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Figure 3: We compare the region favored by PAMELA (green bands) and ATIC (red regions 
within the bands) with the bounds from HESS observations of the Galatic Center [29] (blue 
continuous line), Galactic Ridge [41] (blue dot-dashed), and SgrDwarf [42] (blue dashed) and 
of observations of the Galactic Center at radio-frequencies v = 408 GHz by Davies et al. [51] 
(red lines) and at v ~ 10 14 Hz by VLT [52] (upper purple lines, when present, for equipartition 
and constant magnetic field). We considered DM annihilations into e + e~ (left column), fi + fi~ 
(middle), t + t~ (right), unity boost and Sommerf eld factors and the NFW (upper row), Einasto 
(middle), isothermal (lower) MW DM density profiles and the NFW (upper), large core (middle 
and lower) Sgr dSph DM density profiles. 
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Figure 4: As in the previous fig. 3, but for the cases of DM annihilations into W + W (left), 
bb (middle), ti (right). 
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and radio observations are violated by one or two orders of magnitude. As we conservatively 
fitted the various data-sets, our results are robust. 

Plots are made assuming unit boost factors and Sommerfeld enhancement, and our results 
remain unchanged (up to a rescaling) if these factors are constant. In line of principle, it is 
possible that the e excess (if dominantly due to DM annihilations in the solar neighborhood) 
is significantly more enhanced than DM annihilations around the Galactic Center. However, 
numerical simulations suggest that boost factors much larger than unity are unlikely [16], and 
we have verified that a realistic subhalo population among those discussed in Ref. [55] actually 
produces 0(1) boost factors for ~ 100 GeV positrons and gamma-rays in a 10 -5 cone towards 
the galactic center. It should however be mentioned that a boost factor of a few units is possible 
for positrons, and the first effect of adding substructures to the the smooth DM component 
of the Milky Way actually is to reduce the flux from the Galactic center by a factor (1 — f) 2 , 
where / is the fraction of the mass of the smooth halo that goes into clumps, which is expected 
to be (9(0. 1). In absence of a precise prescription, we limit ourselves to caution the reader that 
this introduces an (9(1-10) uncertainty on the exclusion plots discussed above. Similarly, the 
variation in the Sommerfeld enhancement due to the different DM velocity dispersion at the 
Galactic center can only lead to 0(1) uncertainties on the constraints, since the annihilation 
signal does not arise from regions near the Galactic Black Hole horizon where DM has a larger 
velocity dispersion. 

In order to perform a model-independent analysis, we considered DM annihilations into 
pairs of SM particles. Recently, it was proposed that DM might instead annihilate into some 
new light particle with mass m < m p that can only decay into SM leptons i or pions in view 
of kinematical constraints [22]. Such models lead to a 7-ray flux reduced by about a factor 
\n(M/mi)/hi(rn/rni) ~ 2, only mildly alleviating the HESS bounds considered in this paper. 

The DM annihilation interpretation of the PAMELA/ATIC excesses is compatible with 
gamma and radio observations if the DM density profile is significantly less steep than the 
Einasto and NFW profiles. We also notice that this would mean that e observed by PAMELA 
and ATIC dominantly come from regions of the galaxy that are close to us (as the isothermal 
profile predicts less concentration of DM at the galactic center) and thereby suffer little energy 
losses. This implies in particular that the direct annihilation channel DM DM — > e + e~ is 
disfavored: in this case the spectrum would remain close to a peak at E = M, which seems to 
be disfavored by the PAMELA+ATIC e ± spectra (that show a broader shape). In other words, 
a combined fit of PAMELA, ATIC and photon data from a e + e~ primary channel does not 
yield a good fit for smooth DM density profiles. Other channels that produce broader spectra 
(such as /U + /U~) do a better job. 

Our results are compatible with other recent analyses [54], including those [38] performed 
assuming the very steep DM density profiles of Ref. [28] . 
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